Generation and Eight-port Homodyne Characterization of Time-bin Qubits 
for Continuous-variable Quantum Information Processing 
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We experimentally generate arbitrary time-bin qubits using continuous- wave light. The advantage 
unique to our qubit is its compatibility with deterministic continuous- variable quantum information 
processing. This compatibility comes from its optical coherence with continuous waves, well-defined 
spatio-temporal mode, and frequency spectrum within the operational bandwidth of the current 
continuous-variable technology. We also demonstrate an efficient scheme to characterize time-bin 
qubits via eight-port homodyne measurement. This enables the complete characterization of the 
qubits as two-mode states, as well as a flexible analysis equivalent to the conventional scheme based 
on a Mach-Zehnder interferometer and photon-detection. 
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There are two complementary approaches in optical 
quantum information processing: discrete- variables (DV) 
and continuous- variables (CV). DV experiments are con- 
ducted by qubits represented by single-photon optical 
pulses [lH6|. However, due to inefficient generation and 
imperfect detection of qubits, most of the experiments 
are probabilistic and require post-selection @, HI • In con- 
trast, CV experiments rely on the wave nature of light. 
They can be performed deterministically via quadrature 
squeezing, highly-efficient homodyne detection and feed- 
forward operations, at the expense of relatively low oper- 
ation fidelities 0, [To[ . Recently, there has emerged a "hy- 
brid" approach to combine both techniques to circumvent 
the current limitations [11[. Its major advantage is deter- 
ministic operation of qubits with CV techniques; one of 
the most striking examples is deterministic quantum tele- 
portation of qubits with a CV teleporter, as is proposed 
in Refs. [l2|,[l3|. The recent experiment on CV teleporta- 
tion of highly non-classical optical pulses [14[ opens the 
way to this hybrid teleportation. However, typical qubits 
are generated by pulse-pumped spontaneous parametric 
down-conversion (SPDC) [Ml and thus have no optical 
coherence with continuous-waves on which the CV tele- 
porter is based. Furthermore, the bandwidth of these 
qubits is orders of magnitude wider than the operational 
bandwidth of the CV teleporter (only around 10 MHz). 

Here we overcome this incompatibility by generating 
a time-bin qubit using CW light. This qubit consists 
of two temporally-separated optical pulses, described as 
a superposition of a photon in one pulse |1,0) and the 



other pulse |0, 1) 



c |l,0) + cie**|0,l). Thus 



far, such qubits have been prepared by pulsed lasers 
[H[ and used as a key resource for various DV exper- 
iments over long distances (e.g., quantum cryptogra- 



phy [15|, ll6| , quantum teleportation [17j , and entangle- 
ment swapping [HI). In contrast, our time-bin qubit 
is generated from a CW-pumped nondegenerate opti- 
cal parametric oscillator (NOPO), which is a cavity- 



enhanced version of the SPDC. The NOPO cavity en- 
hances the SPDC process only inside its resonant mode, 
thereby generating qubits in a well-defined and controlled 
spatio-temporal mode. The 6.2 MHz bandwidth of the 
resultant qubit is within the bandwidth of the squeezing 
and homodyne-based feedforward operations (currently 
up to tens of MHz), which mainly limits the bandwidth 
of the CV teleporter in Ref. [l4j|. Thus, our qubit is fully 
compatible with the teleporter as well as more advanced 
CV experiments based on such techniques. Another ad- 
vantage is that our qubits can be readily teleported by a 
single CV teleporter (such as Ref. since two pulses 
have the same polarization, while polarization qubits re- 
quire two C V teleporters (one for each polarization) [l3| . 
Such high compatibility should open the way for further 
hybrid protocols. 

For the characterization of time-bin qubits, we develop 
a scheme via eight-port homodyne measurement 
l2l| . Unlike the other schemes based on photon detec- 
tion Q or regular homodyne measurement [H, 0], our 
scheme provides a simple and efficient way to indepen- 
dently measure two pulses without varying any optical 
phase of the system. This enables the reconstruction of 
the complete two-mode density matrices, not only in the 
qubit subspace spanned by {|1, 0), |0, 1)} but in much 
wider photon-number space including vacuum and multi- 
photon components. Additionally, we show that the 
analysis equivalent to the conventional photon-detection 
scheme ^ can be performed more flexibly for various 
detection conditions with the same setting. 

Our scheme [Fig. DJa)] is based on the setup for gener- 
ating single photons in Ref. [22|. In this setup, an NOPO 
is weakly pumped to produce correlated photon pairs in 



signal and idler modes written as [23|, |2 



dtdt'C(t-t')ai(t)al(t') 10)^. (1) 
Here, a|(t) and al(i) are creation operators in the sig- 
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FIG. 1: (Color online) (a) Schematic, (b) Experimental setup. 
AOM, acousto-optic modulator; APD, avalanche photodiode; 
BS, beam splitter; CW, continuous wave; LO, local oscillator; 
NOPO, nondegenerate optical parametric oscillator; SHG, 
second harmonic generation; and VBS, variable beam splitter. 



nal and idler mode respectively, and C(t — t') denotes 
the correlation function between these two modes. A 
photon-detection event in the idler mode at time t = 
projects this state onto a\(t = 0) 1 0) i . It produces a sin- 
gle photon state A\ |0) s , where A\ = f dtC (t)d\(t) / N 
and N = (J dt \C(t)\ 2 f^. The temporal mode of this 
photon is defined by the mode function f(t) = C(t)/N, 
which converges to fo(t) = y^e -7 ^ in the weak pump- 
ing limit (7: NOPO bandwidth) [HI . 

To obtain time-bin qubits, we introduce an unbal- 
anced Mach-Zehnder interferometer (MZI) into the idler 
channel as in Ref. [5|. This MZI manipulates the time- 
correlation between the signal and idler modes before 
photon detection, thereby reshaping the temporal mode 
in which photons are generated. This scheme originates 
from the theoretical proposal in Ref. [26| and the subse- 
quent experiments in Refs. 27]-29|, m which an unbal- 
anced MZI is introduced to violate Bell's inequality by 
generating energy-time entanglement. We use this en- 
tanglement as a resource of time-bin qubits by detecting 
one of the entangled modes (the idler mode). This MZI 
reduces the generation rate by discarding one of the two 
output ports (it can be avoided in theory as mentioned 
below), but does not degrade the purity of the generated 
qubits. 



To model this scheme, we first introduce an auxiliary 
vacuum mode "v" , and define the beam-splitter (BS) op- 
eration on mode "i" and "v" as replacements of dj(t) —> 
rdl(t)-pe-^at(t) and a+(t) pe^ aj(t)+rat (t), where 
r, p and $ are the amplitude transmissivity, reflectivity 



and relative phase, respectively (r 



(ti, pi, $1 = 0) transforms Eq. (pp) into 



/ 



dtdt' [ Tl C(t-l/)ai(t)a[tf) 

- Pl C(t-t')ai(t)al(t')]\0) s \0) i \0) v . (2) 



By introducing a time delay At in mode "v" , the C(t — t') 
on al(t)al(t') is replaced by C(t - (t' - At)). Mode "i" 
and V are then recombined by BS-2 with (r2,p2>3>2)- 
After this operation, the only term relevant to the photon 
detection on mode "i" has the form 



/ 



dtdt' [nT 2 C(t - t') - p lP2 e i ^C(t - t' + At)] 

alWa^'mjOMO},. (3) 



By projecting this state onto d\(t = 0) |0) i and tracing 
out the unused mode "v" , we obtain 



j dt [nT 2 C(t) - PlP2 e^C(t + At)] ot(t) |0> s . 



(4) 



When At is sufficiently long compared to the photon co- 
herence time (~ I/7) between the signal and idler modes, 
the two modes defined by A\ = j dtC{t)d\{t) /N and 
A\ = fdtC(t + At)d\(t)/N can be regarded as orthogo- 
nal. Thus, from Eq. ((4]), the two- mode state can be de- 
scribed as |^)i2 a r i r 2 |0)2 — Pip2^ 2 |0) x |1) 2 . The 
delay (t f —> t' — At) in the idler mode induces earlier 
time-bin (£—>•£ + At) in the signal mode, and the two 
time-bins constitute a single time-bin qubit. The coeffi- 
cients of the qubit are determined by the splitting ratio 
of both BS-1 and BS-2, as well as the recombining phase 
at BS-2; we experimentally show that they are arbitrarily 
tunable. In the weak pumping regime, the mode func- 
tion of each time-bin is written as fi(t) = y^e -7 ^ and 

/2CO = v /7 e_7|t+A ^ 1 - Note that > when the Photon de- 
tection is performed on the other output port of BS-2 



instead, |?//) 



riP2|l) 1 |0) 2 + p 1 r 2 e^|0) 1 |l)2 ^ ob- 



1). BS-1 



tained. By adding photon-detection events of this port 
and performing an appropriate unitary transformation 
thereafter, the MZI-induced reduction of the generation 
rate is avoidable, though not experimentally shown here. 

Our detailed experimental setup is shown in Fig.QJb). 
Part of the output beam of the Ti: Sapphire laser (wave- 
length: 860 nm) is frequency doubled by a second- 
harmonic generation cavity (frequency: ujo —> 2ujq). Its 
frequency is then shifted with an acousto-optic modula- 
tor (AOM-1) by 590 MHz (2uo 2uo + Ac*;), which cor- 
responds to the free spectrum range (FSR) of our NOPO. 
10 mW of this beam pumps the NOPO [half- width at 
half-maximum (HWHM): 6.2 MHz (= 7/271-)], producing 
signal (coo) and idler (uo+Acj) photons. Two weak coher- 
ent beams (cjo, t^o + Ao;) are also injected into the NOPO 
to lock the subsequent system. Both beams are alterna- 
tively blocked and unblocked by AOM-3 and 4 at a rate 
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FIG. 2: (Color online) Eight-port homodyne analysis of time- 
bin qubits. BS, beam splitter; and MaxLike, maximum like- 
lihood algorithm, (a) Reconstruction of two-mode density 
matrices, (b) The analysis equivalent to the conventional 
photon-detection scheme. 



of 2 kHz. The signal and idler modes are spatially sepa- 
rated by a cavity (HWHM: 22.2 MHz, FSR: 1150 MHz). 
The idler mode is then injected into a MZI after filtering 
out irrelevant NOPO modes via two filter cavities. In the 
MZI, part of the idler mode is split off by a variable BS 
(VBS), composed of a half- wave plate and a polarizing 
BS. The split-off part is transmitted through a 50m opti- 
cal fiber (PM780-HP, Thorlabs), and is then recombined 
with another VBS. Here the delay is set to At = 242 ns 
so that the overlap between two time-bins is negligible 
(/ h(t)f 2 (t)dt ~ 0.08%). The optical path length differ- 
ence is passively stabilized by surrounding the fiber with 
heat insulating material, as well as actively locked to a 
target value with piezo actuators. One output port of 
the MZI is sent to an avalanche photodiode (APD) to 
herald time-bin qubits, whereas the other is monitored 
by a photodetector to phase lock the MZI. 

In general, complete homodyne characterization of 
two-mode states requires two-mode quadrature measure- 
ments at various phase sets (0i,6 2 ) [6]. One possi- 
ble way to carry out this in our case is to spatially 
separate two time-bins and then perform a homodyne- 
measurement on each bin with careful control of the 
local oscillator's phase. However, since two bins have 
the same polarization and are temporally close to each 
other, such a separation requires high-speed switching 
with high precision. In Ref. [5[, only the relative phase 
of two time-bins are scanned instead of fully scanning 
local oscillators' phases, but this scheme is inapplica- 
ble when characterizing unknown qubit states. Here, 
we propose a more efficient scheme via eight-port homo- 
dyne measurement 19|-2l| where no optical phase scan- 
ning is required. First, both time-bins are split with a 
50/50 BS for two homodyne measurements for orthogonal 
quadratures x and p. By multiplying fi(t) and f2(t) to 
each homodyne photocurrent, we can obtain the quadra- 



tures (xi,pi,X2,P2) with extra vacuum noise from the 
unused port of the BS [Fig. [2(a)]; this vacuum noise is 
the inevitable consequence of the simultaneous measure- 
ment of conjugate variables x and p [30(] . From this set, 
quadratures at any phase can be calculated as Xj (0) = 
Xj cos# + pj sin# (j = 1,2). Therefore, from one set, we 
can randomly select a phase set (0i,0 2 ) an d construct 
two independent tomography data sets (x\ (#i), £2 (#2)) 
and (x\{6\ + 7r/2), £2 (#2 + 7r /2)) which refer to orthog- 
onal quadratures. In the experiment, 100000 quadrature 
sets are acquired and 200000 tomography data sets are 
constructed. The maximum likelihood algorithm (3l| is 
then used to reconstruct two-mode density matrices in 
the Fock basis: p = Y^T,i,m,n=o Pkimn \k,l) (m,n|. Here 
we compensate only the extra vacuum noise added in the 
dual- homodyne measurement, and not other experimen- 
tal imperfections. The scheme mentioned above can be 
easily extended to single- or higher-mode states, which 
shows that eight-port homodyne measurement is a use- 
ful tool for characterizing multi-mode states. Note that, 
in the present setting, the eight-port homodyne measure- 
ment is directly measuring the two-mode Q function of a 
qubit [HI, [HI . Once the Q function is estimated from the 
collected data, the corresponding density matrix can be 
calculated in theory Instead, we have mapped the 
measured data (xi,pi, x 2: p 2 ) to (#1, xi(0i), 2 , x 2 (0 2 )) : 
from which we can directly calculate the density matrix 
using the standard maximum likelihood algorithm (3l| . 

We generated eight types of qubits: ( 1 1 , 0) + 
j*\0,l))/y/2 and (2|1,0) + e** |0, l))/>/5 with $ = 
0,7r, ±7r/2. Figure [3(a) shows the quadrature variance 
at each time point calculated from the 100000 quadra- 
ture traces where the vacuum variance is normalized 
to 1/2 (h = 1). Two time-bins appear as the two 
peaks of the trace, and the peak height above the vac- 
uum noise level corresponds to the probability of ex- 
isting photons (23|. Experimental density matrices in 
Fig. [3(b,c) show not only the qubit submatrix spanned 
by {|1, 0) , |0, 1)} but also the vacuum and multiphoton 
contributions. Ideally, generated states have components 
only in the qubit submatrix. In practice, the density ma- 
trices show 18 ± 1% of a vacuum (poooo), 77 ± 1% of a 
qubit (pioio + P0101), and 5 ± 1% of multiphoton compo- 
nents (1 — poooo — P1010 — Poioi)- The ratio of diagonal 
elements pioio to poioi is equal to the ratio of |1,0) to 
|0,1), whereas the off-diagonal elements piooi and pono 
demonstrate the superposition of |1,0) and |0, 1) at the 
target phase. The encoded quantum information can be 
read out by extracting and renormalizing the qubit sub- 
matrix. The average fidelity of each submatrix with its 
target state is 0.989 ± 0.004, showing the highly precise 
qubit preparation. 

The 18% vacuum contribution is well explained by 
the estimated loss of 1 — 77^,11 = 16%, where ry a n = 
^NOPO^vis^pr^det^APD is the overall efficiency. Here 
^nopo = 0-98 is the escape efficiency of the NOPO, 



^vis = 0.98 the average mode- matching visibility at eight- 
port homodyne measurement, r] pr = 0.96 the propaga- 
tion efficiency from the NOPO to the homodyne detec- 
tors, 7/det = 0.95 the average detection efficiency given by 
the quantum efficiency of photodiodes and the electronic 
noise, and t?apd = (Ctot - Cdark)/Ctot = 0.98 the purity of 
the photon-detection event given by the total count rate 
Ctot = 5800 s _1 and the dark count rate Cdark = 80 s _1 . 

Our eight-port homodyne setup also enables the anal- 
ysis equivalent to the conventional detection scheme in 
Fig. HJa), where two time-bins are recombined via a 
MZI and then the existence of photons at each output 
mode is monitored by photon detections Dl and D2. 
In this case, Dl and D2 correspond to the projection 
onto orthogonal qubit bases, and if the BS parameters 
are chosen properly, photons arrive at only one of the 
two. In our system, this BS can be realized by analyti- 
cally mixing measured quadratures of two time-bins after 
the measurement [Fig. [5(b)]. For a quadrature data set 
(xi,pi,X2,P2), the complex amplitude of each mode can 
be written as a, = ( Xj +ip )/V2 (j = 1,2). When 
the BS recombines these two time-bins with (r', //,</>), 
the output amplitude (a[,a 2 ) can be calculated as a[ = 
t'ol\ + p'e l ^a2 and a! 2 = — p'e~ l ^a\ + r'o^. Finally the 
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FIG. 4: (Color online) (a) Conventional photon-detection 
scheme, (b) The <j) dependence of single photon components of 
two modes which are subject to photon detections Dl and D2. 
(c,d): Experimental Wigner functions W(x,p) and photon- 
number distributions P(n) at <j) — tt/2. 




FIG. 3: (Color online) Experimental results. Quadrature 
variance at each time point (a) and density matrices of (b): 
(|1, 0) - i|0, l»/>/2 and (c): (2|1, 0) + |0, 1))/VE. 



output quadratures at any given phase 6 can be cal- 
culated from these amplitudes as x'j(0) = (pi^e~ %d + 
a f *e z0 )/y/2. Therefore, we can compose a data set 
(^1(^1)5^2(^2)) fc> r any desired BS parameters (r', //,(/>) 
and any phase set (#i,#2)- As a result, the density ma- 
trix of each mode subject to Dl or D2 can be recon- 
structed by maximum likelihood algorithm compensating 
the extra vacuum noise. In the special case of choosing 
(t',//,0) = (co,ci,— $) or (ci,co, — 3> + 7r) for an initial 
time-bin qubit \^) 12 = Co |0) 2 + C\e % ® |0) 1 |1) 2 , the 
two modes would theoretically be decomposed into a sin- 
gle photon state and a vacuum state. In our eight-port 
homodyne scheme, the whole set {(r 7 , //, (/>)} of the Dl- 
D2 measurement results can be obtained from a single 
quadrature data set. Some related discussions on this 
photon- vacuum decomposition were made in terms of the 
cross section of two- mode Wigner function in Refs. [H|,[6[. 

We analyzed the data set of the qubit = (|1,0) — 
z|0,1))/a/2 for a fixed r' = p' = l/y/2 and varying <p. 
The <p dependence of the single photon component of Dl 
and D2 [Fig. g^b)] shows the fringe visibility of 96 ± 2%. 
It reaches its minimum and maximum around ±7r/2, 
demonstrating the generation of the target qubit state. 
The reconstructed states at (j) = tt/2 are explicitly de- 
composed into a single photon and a vacuum. Figure |4jc) 
shows P a (l) = 0-79 of a single photon component and a 
strong negativity of Wigner function VK a (0,0) = —0.187. 
Figure H(d) is an almost pure vacuum with Pb(0) = 0.98. 
Pb(l) = 0-02 of a photon still remains, possibly because 
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the BS parameters (t',//,0) are not optimal and the 
phase $ of the qubit fluctuated during the measurement. 

In summary, we have generated arbitrary time-bin 
qubits in a well-defined spatio-temporal mode by using 
a CW-pumped NOPO. The generated qubit is compati- 
ble with the current CV operations, thereby providing an 
essential tool for hybrid approaches in quantum informa- 
tion processing. Our eight-port homodyne scheme effi- 
ciently reconstructed the full two-mode density matrices 
of the qubits, and also enabled the analysis equivalent to 
the conventional photon-detection scheme by a numerical 
BS operation. This characterization scheme is universal, 
and applicable to single- or multi-mode states encoded 
in different temporal modes , spatial modes 0, E[ or 
polarization modes 
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